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ABSTRACT 

A group of stations in the North American Arctic region have been 
analyzed for statistical determination of temperatures at mandatory 
pressure levels p.. For each station the temperature at a key level 
(called the forciAg-level temperature) peculiar to the station has 
been forced in at the first step, and retained at each subsequent step 
in the development of the stepwise regression equations giving tem- 
perature at the mandatory levels. In general, eight-step prediction 
equations in terms of other temperatures in the vertical were found to 
give specification of T(p,), with percentage explained variance of 
close to 0.99. As a result of this definitive property, the best- 
estimate of the regional atmosphere which is conditionally dependent 
upon the existence of an extreme probability of the forcing^level 
temperature is obtained with a high degree of confidence. 
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1 . Introduct ion 



A preliminary set of model atmospheres giving both the warm and cold 
extremes at certain reference elevations with a 1 ^ level of expectancy 
has been proposed as the Preliminary M1L-STD-210B by ETAC (Environmental 
Technical Applications Center of the USAF). In essence this model of 
cold and warm extreme temperatures is based upon listing the coldest 
(warmest)temperatures at each of the eight standard levels listed in 
Table 1. These extreme temperatures are based upon observations of the 
location of extreme temperature noted on global temperature maps at the 
level noted. The extremes are observed to occur at the geographically 
diverse locations listed in Table 1. 

The values of the extremes in Table 1 were computed from the assumed 
existence of a normal distribution of temperature in both the coldest 
(January, N. Hemisphere) and warmest (July, N. Hemisphere) months. Thus 
these extremes have been computed based upon the expectancy (true for a 
normal distribution) of finding 1^ of all temperature-departures from the 
sample mean at each location and level at 2.3267 times the standard 
deviation O at that level and location. In the process of making the 
estimation of temperature- extremes , no allowance has been made for the 
possibility of inter-level correlations which have been shown to exist 
at typical stations (Cole and Nee, 1965). 

Because of the neglect of vertical consistency in assembling the tem- 
perature-extreme data for the Preliminary MIL-STD-210B climatic model, an 
unrealistic range of temperatures is required in design of equipment 
needed to operate in these atmospheres. As a result, it has been proposed 
that the geographical sites which are proposed in Table 1, as correspond- 
ing to world-wide temperature extremes at the given levels^ be used to 
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TABLE 1. Locations of proposed extreme temperatures in the 
cold world-wide and warm world-wide cases (after ETAC). 





Cold-Extreme 


Warm- Extreme 


Level 


Location 


January 

Mean 

(°C) 


STD. 

DEV. 

(°C) 


Extreme 

(°C) 


Location 


July 

Mean 

(°C) 


STD. 

DEV. 

(°C) 


li 

Extreme 

(°C) 


SFC 


Oymyakon, 

USSR 


-50.6 


4.4 


-60.6 


Insalah, 

Algeria 


40.9 


3.9 


49.0 


850mb 


Oymyakon, 

USSR 


-35.1 


5.5 


-hl.Q 


Insalah, 

Algeria 


28.9 


2.2 


34,0 


700mb 


Hall Beach, 
NWT 


-27.3 


6.5 


-42.4 


Babylon, 

Iraq 


17.0 


2.2 


22.1 


500mb 


Resolute, 

NWT 


-43.0 


4.6 


-53.2 


New Delhi, 
India 


-4.3 


3.9 


4,0 


300mb 


Thule, 

Greenland 


-60.6 


2.2 


-66.0 


New Delhi, 
India 


-25.8 


3.9 


-16.0 


200mb 


Thule, 

Greenland 


-59.7 


7.2 


-75.0 


Alert , 
NWT 


-42.3 


2.5 


-36.5 


150mb 


Karachi, 

Pakistan 


-65.0 


6.0 


-79.0 


Alert, 

NWT 


-43.3 


2.5 


-37.5 


lOOmb 


Singapore 

1 


-83.0 


3.0 


f 

-90.0 

L 


Thule, 

Greenland 


-43.8 


2.2 


-37.2 
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provide ’’forcing temperatures" at the level. Based upon inter-level tem- 
perature-statistics at each station, the most likely temperatures at the 
eight additional levels (the 70 mb level is to be included) in the atmos- 
phere are then to be estimated with optimal accuracy. 

It is to be noted that Table 2 lists only six different stations in 
definition of the cold-extreme, and five in definition of the warm-extreme 
atmosphere. Hence, a real-time data analysis would require a detailed 
climaticalogical examination of 11 different stations for the purpose of 
establishing regional extremes. 

2. A statistical model for regional extreme atmospheres. 

As indicated in the study by Cole and Nee (1965), high values of 
inter-level simple correlations of temperature frequently exist. It was 
decided in this study to obtain an initial data file on magnetic tape for 
the two-a-day rawinsonde stations (all in the North American sector) 
identified by code symbols Cl and Wl in Table 1. These symbols are 
employed to suggest the cold and warm atmosphere feasibility study. The 
data was provided by the Environmental Data Service, National Climatic 
Center, through the financial assistance and kind cooperation of the Com- 
mander, Naval Weather Service. 

The other stations in Table 1 not bearing the code symbols Cl or Wl 
have also been provided through the data sources just listed. However, 
these secondary data records are not available in the convenient summary 
form, nor in such abundant population samples as those stations which were 
coded Cl or Wl so that more care is required in the data-handling of these 
secondary stations. This second set of data-stations is the subject of 
ongoing work which will appear as a later report in this series. 
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For the stations marked code symbols Cl and Wl, a stepwise multiple 
regression technique was employed. This program, known as BIMED 02R, is 
available in the Program Library of the W. R. Church Computer Center of 
the Naval Postgraduate School. The temperature data was arranged in a 
decreasing pressure sequence but with standard pressure-spacings as fol- 
lows (Table 2): 

TABLE 2, 





J=1 


J=2 


J=3 


J=4 


J=5 


J=6 






lOOOmb 


950 


900 


850 


800 


750 




II 


700 


650 


600 


550 


500 


450 


J = 12 


J = 13 


400 


350 


300 


250 


200 


175 


J = 18 


J = 19 


150 


125 


100 


80 


70 




J = 23 



Actually, each Cl and Wl radiosonde was listed with both a height and 
temperature at each of these levels. Over the four year data sample 
(1967-70), these stations provided the following population samples: 



C-1 Stations 

Hall Beach (700mb) 243 samples 
Resolute (SOOmb) 237 samples 

Thule (300mb) 213 samples 

Thule ....•(200mb) 213 samples 



W-1 Stations 

Alert (200mb) 244 samples 

Alert •..••(ISOmb) 244 samples 
Thule (ISOmb) 236 samples 



The program BIMED 02R permits for each station, the specification of 
a different dependent variable or predictand in each subproblem. Also, 
the number of independent variables to be tested for admission to regres- 
sion equation can be limited. In the problems considered here an upper 
limit of eight independent variables was found to be given optimal speci- 
fication of each of the dependent variables. Another convenient feature 
of the BIMED 02R program is that the forced entry into the regression of 
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the temperature value being tested as an extreme (the so-called ” forcing 
temperature’*) is allowed, even though its contribution to the explained 
variance may become small. For example, in the case of Resolute, which 
corresponds to a winter extreme at 500 mb, the multiple regression equa- 
tions at the various ’’mandatory levels” were set up in the forms: 



\ - A + B + C (1) 

Here is the 8 -level predictand chosen alternately as 

rp . •T' rn rp rp rp rp rp rp 

1* 4* 7* ^15* 17* 19* 21* 23 

where the subscript M is the temperature identifier associated with the 
pressure level indicated in Table 2. In summary, the predictands 
desired are 

the temperature at 1000 mb 

T^, the temperature at 850 mb 

Ty, the temperature at 700 mb 

^ 11 * temperature at 500 mb 

the temperature at 300 mb 

the temperature at 200 mb 

the temperature at 150 mb 

^ 21 * temperature at 100 mb 

T 23 * the temperature at 70 mb 

The form represents here a linear combination of up to 7 independ- 
ent variables, apart from the forced temperature in the case of 
Resolute, NWT. The variable varies from case to case, as the predictand 
is changed from T^ to T^, to Ty, etc. X^ usually involves temperature 
predictors which are physically close to the level being specified as 
the dependent variable T^^. For example, at Resolute, T^ [= T(850)] 
was best- fit by the following expression: 
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T 



4 



3745 - .0849 T 
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+ 1.0 X 



% 



X, = .3614T_ + .85771^ - .2892T, + .0819T, 
4 3 5 6 7 



+ .0902T„ + .0598T,,- .0486T,_ 
9 14 17 



with a multiple regression coefficient R = .9861. In general, somewhat 
similar results applied at each of the other mandatory pressure levels 
1000, 700, 500, 300, 200, 150, 100, 70 mb except that the multiple R 
expressing Tj^[= T(IOOO)] and T(300)] in terms of and X^ were 

somewhat lower than the remaining correlations in the set. Similar 
statements applied to all stations at both extreme-times January and 
July. 

The function X^^ is called in this study a "structure-function" in 

that to a large degree it accounts for most of the variation of the T 

M 

being predicted, even though the part played by the forcing variable 
is by no means negligible, 

A least squares fit of form Eq* (1) may be shown to produce residuals 
(computed against the observed values) which are distributed normally 
relative to the regress ion- value, together with a variance given by the 



variable regressions tested here, most of the multiple regression coeffi- 
cients had values well in excess of ,99 (with the exception of the two 
levels already noted, 1000 mb and 300 mb), so that a near functional 
relationship existed between the predictor and the predictand set. 

The specification statistics for the different stations in the North 
American sector of Table 1 are listed in Section 3, 



standard 




In the set of dependent 
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3. Regression-specification of mandatory- level temperatures using the 
forcing- level temperature. 

The stepwise regression procedure of BIMED 02R was employed with the 
appropriate forcing-level temperature required at the first step of the 
regression and retained in the selection of the following seven predictors. 
These seven predictors were allowed to enter the stepwise regression in 
accordance with the programmed requirement that at each step k, the variable 
added explained a maximum of the unexplained variance remaining after (k-1) 
selections (k=2, 8) . 

Tables 3a, 4a,..., 9a show the results of the stepwise regression applied 
to the specification of the temperatures T^, T^, T^, 

^21’ ^23 mandatory levels by an equation of form (1), using an eight- 
predictor equation for each T^ just listed. Table 3a lists for the January 
fulj. data-sample at Hall Beach the mean, standard deviation, multiple cor- 
relation coefficient and standard error at each level other than at 
forcing level. Table 3b lists analogous results for a nominal 10 ^ of the 
January data sample corresponding to the 10^ cold extreme sample. 

Table 4a and 4b list similar results for Resolute in January with 
T^^[= T(500)] as the forced variable. Tables 5 and 6 list analogous 

results for Thule, using T^^ and T^^ as forcing variables, respectively. 

2 

It is to be noted that the fractional explained variance R of the predic- 
tand in each case, does not fall of appreciably, in general, as one pro- 
ceeds along any row of these tables from part (a) to part (b) . 

Tables 7, 8, 9 similarly list the July-warm statistics, using respec- 
tively the full data-samples at Alert [with forcing temperature T^= T(200)], 
Alert [T = T(150)] and Thule [T = T(IOO)] in the left half of the tables. 

J J 

The prediction results corresponding to the nominal 10^ of the warmest 
July forcing-level temperatures appear on the right half of these tables. 
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The selection of nominal 10^ extreme cases was made by separating out of 
the full sample all cases where the forcing temperature T^. lay in the ranges 
T < T - 1.2817a for cold extremes 

J vl 

Tj > Tj + l*2817aj for warm extremes 

The criteria used in establishing these 90^^ nominal extremes was that the 
full data samples were normally distributed, a valid first approximation. 

The actual number of extreme data cases is listed at the top of Tables 3b, 
9b, for each case under consideration. 

In order to determine the multiple correlation coefficients and standard 
errors for the right-side 10^ extreme samples of Tables 3b,..., 9b, a regres- 
sion equation of identical form to that developed for the full data-base 
predictands was generated for testing the 10^ extreme-case data samples. 

The test-predictors for the extreme-data were then considered to be the 
forcing variable T^ and the structure-variable X of form determined for 
the full-sample cases. The high values of the multiple correlations thus 
determined for the January extreme-data cases serves to verify the concept 
that the predictands corresponding to the extreme-data may be anticipated 
through a joint knowledge of the forcing variable T^. and of the structure 
function. The somewhat smaller correlation coefficients for the July data 
is due to the much smaller values of variance in summer, particularly in 
the lower stratosphere. 

A listing of the nominal extremes of the forcing temperature for other 
probability levels of interest, e.g., the 5^, the 1^ was also made assuming 
the existence of a Gaussian distribution of temperatur^ and appears in 
Table 10. A separation of the 10^ nominal extremes of T^ was then made 
wherever possible into the classes suggested in Table 10, that is beyond 
the .01 probability level; in the probability range .05 to .01 and in the 
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range .10 to .05. The results of these probability stratifications on the 
predictands at T^, T^, T^, ^17’ ^^ 21 * ^^23 discussed in Section 4. 

These additional critical extremes are given by 



Tj(.05) = Tj ± 


1.640 a 


(2) 


Tj(.Ol) = T j + 


2.3267a 

J 


(3) 



Equations 2 and 3 tend to give excessive estirmtes of the critical 
extremes when the observed Q is too large to be representative of a normal 
distribution. This happens primarily in the winter season at low strato- 
spheric levels, e.g., Thule at 300 mb. At this locale there tends to be 
bimodal alternations between warm and cold stratospheric regimes and an 
unrepresentatively large ^ results from this sample. The ETAC estimate 

sj 

(Table 1) of the extreme temperature at Thule at 300 mb in winter also 
seems to reflect this same difficulty. 

The following temperature classes were then considered as input data at 
the Cl and Wl station forcing levels, respectively 

Cold classes: (a) (b> T T (c) T T jq 

Warm claasaa , (a) T (b) T_„j>^>T ,,5; (c) T„5>T>T 

In most cases studied, reasonable sample-sizes existed in each class (a), 
(b), (c). However, in several cases, the extremes did not reach the 1^ 
nominal-value predicted from an assumed Gaussian distribution. 
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TABLE 3. Regression statistics at mandatory pressure levels at Hall Beach using 
Tj= T(700) as the forcing-level temperature. Part (a) refers to the full-data 
January sample; part (b) refers to the nominal 10^5 cold extreme sample of T(700). 





M-l ^ 
O O 




















O 

u ^ 
o 

u <u 


















U) 

<D 

W 

cd 


U 4 J 
w Cd 


I— t 


CO 




CTN 


<t 


00 


1—1 


r-- 


G 


0 


uo 




UO 


as 


uo 


UO 


00 


• 

X> 4 J 


C 30 

• 






vD 

• 


0 

• 


00 

• 


1—1 

• 


00 


o 


4 J CO 

W W 


CN 








t— 1 




1—1 




vD 




















CSJ 




















II 


CU 




















i-H • 

D.I »~l • 


00 






0 


00 


00 




1—1 




•H (U 4-1 


CO 


LO 




UO 




r^ 


s£> 


vt 




4J ^ 4-1 


CM 






r-- 


00 


00 


00 


<Js 




1— 1 (U 




(j\ 




CJN 


(J^ 


C 7 ^ 


as 


as 




3 0 0 
:s C-) cj 






) 


• 


• 


• 


• 


• 




• 

> 




















(U 


CO 


m 


vD 




uo 


as 


vD 


r- 




T 3 


G\ 




00 




vD 




CJ^ 






• C-> 


00 


CM 


j— 1 


00 


<JS 


CM 




00 




TJO 


• 


• 


• 


• 


• 


• 


• 






4 J 


CO 


CO 


CO 


CM 




UO 


vO 






CO 




0 


00 


CO 


UO 


00 


0 


1—1 






0 


vO 


sD 


C 7 ^ 




00 




1 — 1 






• 


• 


• 


• 


• 


• 


• 


• 




c 3 CJ 


uo 


m 


as 


UO 


uo 


UO 




1—4 




(DO 


CO 


CO 


CO 










UO 




X 


1 


1 


1 


1 


1 


1 


1 


1 




»— 1 


G 


















(D 


0 


0 


0 


0 


0 


0 


0 


0 




> 


0 


uo 


0 


0 


0 


0 


uo 


0 




CU 


0 


00 


r-- 


uo 


CO 


CM 


1—1 


1—4 






»— 1 


















4-1 

O CJ 




















o 

u ^ 
o 

U CU 
)-4 4-1 


CO 






MD 


1—1 


CO 


VO 


UO 


U) 


w cd 


00 


CM 




as 


0 


00 


CO 


CO 


Q) 


B 


CO 


00 






CM 


00 


00 


00 


U) 

cd 


* ‘l-l 
4 J 


• 

CO 






• 


# 

i-H 






• 


o 


4 J to 
CO W 


















CO 








































Cvj 




















II 


CU 

• 


















13 


CI4 r- 1 • 

CU 44 

I 1 LH 


CM 






CM 


CO 




0 


CO 


x'-N 


(J^ 


CO 




uo 




as 


CM 




cd 


r-i U (D 

500 

sou 


o^ 


C 7 ^ 




as 




00 


<J^ 


<JS 




00 


as 




<JS 


CJ^ 

• 


C 7 ^ 


<ys 


<JS 




• 

> 




















CU 

FtH 


CM 


00 




<r 


00 


1—1 




as 




W 

• 0 

x>o 


0 


CO 


vD 


<ys 


r>. 




1—1 


00 




vO 


CM 

• 


• 


as 

• 


uo 

• 


0 

• 


uo 


vO 




4 J 

CO 






vO 




UO 


vD 


vO 








l— 1 


CO 


m 


CO 


as 


CO 


1—1 


C 7 ^ 






1— I 


vO 


CO 


as 


0 


uo 




<t 




M 




• 


• 


• 


• 


• 


• 


• 




Cd CL? 


vD 


CM 




as 


uo 


CO 


CO 


UO 




cuo 

s 


CM 

1 


CM 

1 


CM 

1 


CO 

1 


uo 


uo 

1 


uo 


uo 




T— 1 


•i 


















CU 


0 


0 


0 


0 


0 


0 


0 


0 




> 


0 


UO 


0 


0 


0 


0 


UO 


0 




CU 


0 


00 




uo 


CO 


CM 


1—4 


1—4 




hJ 


t— 1 

















12 



-57.38 8.942 .9946 .944 70 -54.41 8.936 .9936 1.048 



TABLE 4. Regression statistics at mandatory pressure levels at Resolute using 
Tj= T(500) as the forcing-level temperature. Part (a) refers to the full-data 
January sample; part (b) refers to the nominal 10 % cold extreme sample of T(500). 
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TABLE 5. Regression statistics at mandatory pressure levels at Thule using 

T(300) as forcing-level January temperature. Part (a) refers to the full- 
data sample; part (b) refers to the nominal 10^ cold extreme sample of T(300). 
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TABLE 6. Regression statistics at mandatory pressure levels at Thule using 
Tj= T(200) as forcing-level temperature. Part (a) refers to the full-data 
January sample; part (b) refers to the nominal 10^ cold extreme sample of T(200). 
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TABLE 7. Regression statistics at mandatory pressure levels at Alert using T = T(200) 
as forcing-level temperature. Part (a) refers to the full-data July sample; “^part (b) 
refers to the nominal 10^ warm extreme of T(200). 
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TABLE 8. Regression statistics at mandatory pressure levels at Alert using T = T(150) 
as the forcing-level temperature. Part (a) refers to the full-data July sampfe; part 
(b) refers to the nominal 10^ warm extreme sample of T(150). 
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TABLE 9. Regression statistics at mandatory pressure levels at Thule using T = T(IOO) 
as the forcing-level temperature. Part (a) refers to the full-data July sampfe; part 
(b) refers to the nominal 10^ warm extreme sample of T(IOO). 
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TABLE 10. Nominal temperature-extremes at the indicated 
stations and pressure levels at the 10^, 5^ and 1^ extreme 
values according to a Gaussian distribution. 



Station 


T 

.10 


T 

.05 


T 01 (°C) 


Forcing Level T^ 


90^ extreme 
sample size 


Hall Beach 


-35.6 


-37.9 


-42.3 


L 


= 


T(700) 


26 


Resolute 


-47.1 


-48.5 


-51.3 




= 


T(500) 


25 


Thule, 


Wi 


-62.0 


-63.8 


-67.3 


L 


= 


T(300) 


13 


Thule, 


Wi 


-63.5 


-66.0 


-70.6 




= 


T(200) 


18 


Thule, 


Su 


-41.3 


-40.8 


-39.8 




= 


T(IOO) 


18 


Alert, 


Su 


-39.7 


-39.1 


-37.9 


"j 


= 


T(150) 


22 


Alert, 


Su 


-38.6 


-37.6 


-35.7 




= 


T(200) 


18 
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4. Temperature-estimates at mandatory levels corresponding to 1^ extremes 
of forcing-level temperatures. 

a. Winter extreme atmospheres 

Corresponding to the four winter extreme cases, the multiple regres- 
sion methods which led to Tables 3, 4, 5, 6 also lead to the results of 
Tables 11, 12, 13, 14. In the latter tables, the nominal 10^ extreme sample 
is decomposed into the subsample extreme-classes (a), (b) , (c). These 
classes were defined in Section 3, but their definition is also implicitly 
given in the top row of each table (e.g., Table 11: (a) T (b) T - 
T (c) ^ 05 ~ ^ 10 * Corresponding to the stratified data-sample at each 
forcing level there exist also conditional data-sets at each of the manda- 
tory levels above and below the forcing level. In virtually all cases, 
the mandatory- level data sets in classes (a), (b), (c) exhibit well-defined 
trends between the class-mean temperatures at all levels. This discovery 
afforded credibility to the listed extreme values of the class (a) mandatory- 
level results of Tables 11, 12, 14. Note that in each stratification (a), 
(b), (c), the same sample size exists at each level within each class re- 
gardless of the level under consideration. 

As noted in Section 3, no nominal 1^ data class was realized in 

Table 13; however, the class (b) sample T^. = T(300) comprised a de facto 

extreme 1^ data-sample by actual count. In the comparison of T, with T at 

D c 

mandatory levels over Thule (Table 13), there existed consistent 
inter-class temperature differences, again at all levels. In all cases, 
the reason for these consistent inter-class temperature deviations is that 
the mandatory level temperatures have been found to be strongly correlated 
with the forcing-levels temperature through the multiple regression 
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[ Level 

1 (mb) 

1 

1 


Temperature -means (°C) 
in classes 


Std, deviations 
relative to class 


(°C) 

-mean 




Ca) 
T ^ 


(b) 
T „,-T 


(c) 
T „^-T 


a 


a 


0 




.01 


.01 .05 


.05 .10 


a 


b 


c 


1000 


-35.24 


-35.22 


-34.62 


2.588 


4.047 


3.876 


850 


-39.52 


-35.71 


-32.34 


1.421 


2.475 


i.943 


* 

700 


-44.96 


-39.55 


-36.60 


1.746 


1.175 


0.799 


500 


-44 . 04 


-45.90 


-47.32 


0.948 


3.363 


1.233 


300 


-41.02 


-45.95 


-47.41 


0.937 


5.472 


3.315 


200 


-42.96 


-46.22 


-47.16 


0.924 


6.064 


4.327 


150 


-45.20 


-47.32 


-48.91 


0.767 


7.804 


6.220 


100 


-48.64 


-51.40 


-52.25 


1.284 


8.999 


6.849 


70 


-50.28 


-55.47 


-55.42 


2.047 


10.185 


6.527 


Number 

in 

class 


5 


13 


8 









TABLE 11. Temperature means at mandatory levels at Hall Beach (winter), 
corresponding to temperature extremes of ^700) at the (a) 1^ class, 

("b) 1-5^ class and (c) 5- 10^6 class of probability. The standard devia- 
tions 0 , a , a of the observed temperatures within each class are shown 
in the ® ° right half of the table. 
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Temp.rC) 



FIG. 1. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Hall Beach, NWT, corresponding 
to the class (a) set of 1^ cold extreme occurrences of T-(700) of 
Table 11. The thin solid line depicts the January mean Hall Beach 
vertical sounding (1967-70). The dashed curves the world-wide Pre- 
liminary MIL-STD 210 extreme atmospheres. 
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Level 

(mb) 


Temperature-means (°C) 
in classes 


Std. deviations 
relative to class 


(°C) 

-mean 




(a) 

T 


(b) 

T ...-T .. 


(c) 

T -T . .. 


a 


a, 


0 




.01 


.01 .05 


.05 .10 


a 


b 


C 


1000 


-39.00 


-38.77 


-33.31 


0.748 


2.125 


3.718 


850 


-35.50 


-31.57 


-30.14 


0.698 


1.463 


3.038 


700 


-40.37 


-36.43 


-34.45 


0.742 


1.308 


2.274 


* 

500 


-51.73 


-49.70 


-47.59 


0.173 


0.782 


0.342 


300 


-54.03 


-56.20 


-60.48 


1.225 


2.027 


4.070 


200 


-55.20 


-55 . 24 


-59.15 


0.648 


2.280 


3.561 


150 


-59.37 


-58.09 


-61.12 


0.881 


2.489 


3.098 


100 


-65.70 


-61.00 


-64.59 


0.860 


3.745 


3.248 


70 


-70.77 


-62.50 


-65.88 


0.480 


5.748 


4.064 


Number 

in 

class 


3 


7 


15 









TABLE 12. Temperature means at mandatory levels at Resolute (winter), 
corresponding to temperature extremes of T (500) at the (a) 1^ class, 
(b) 1-5^ class and (cl 5-lO'jl) class of probability. The standard devia- 
tions a , O, , O of the observed temperatures within each class are 
shown ^ ^ in the right half of the table. 
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FIG. 2. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Resolute, NWT, corresponding to 
the class (a) set of 1^ cold extreme occurrences of Tj(500) of 
Table 12. The thin solid line depicts the January mean Resolute 
soimdlng (1967-70) . The dashed curves are Identical to those of 
Figure 1. 
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relationships summarized in Tables 3 through 9. These relationships also 
apply across classes (a), (b), (c), and give rise to the inter-class dif- 
ferences found in Tables 11,..., 17. 

The listings in Tables 11, 12, 13, 14 of the class (a) set of 
temperatures [or the (b) -class if class (a) is non-existen^ at mandatory 
levels makes possible the presentation of the results in graphical form. 

Thus the expected (cold) temperature data corresponding to the 1^ forcing 
temperature appropriate to Tables 11, 12, 13 and 14 have been reproduced 
as Figs. 1, 2, 3, 4. Here the heavy solid line depicts the cold or 1/^ 
”extreme**atmosphere in each case, whereas the thin solid line depicts the 
January (1967-70) mean sounding in each case. Note that the temperature 
deviations between *’mean” and ‘‘extreme” in Figs. 1, 2, 3, 4 are indicative also 
of the class (a) to (c) deviations at the mandatory levels of Tables 11, 



12, 13, 14. 

b. Warm extreme atmospheric cases. 

The procedure in these cases consists of an analysis of the nominal 
10^ extreme-warm July cases at Alert = T(200)], [T^ ~ T(150)]^ and at 

Thule [Tj = T(IOO)] after these extremes have been subdivided into the 
nominal probability classes (a), (b), (c), similar to the procedure of 
Section 4(a) . 

In two of the three cases, the class (a) nominal probability ’^< 

T did not occur (see Tables 15 and 17), but an actual case count reveals 
. u i. 

that in either case there were 2 to 4 data-samples in class (b) . These 
samples were taken as representative of the actual 1^ extreme warm situa- 
tion at the respective forcing levels, and also at the regression- 
dependent mandatory levels. 
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Level 


Temperature-means (^C) 


Std. 


deviations 


(°C) 


(mb) 


in classes 


relative to class-mean 






(b) 


(c) 












T -T 

.01 .05 


T -T 

.05 .10 




^b 


a 

c 


1000 




-25.20 


-20.27 




2.371 


4.942 


850 




-24.83 


-21.87 




1.637 


3.736 


700 




-29.73 


-28.36 




5.287 


3.858 


500 




-42.90 


-42.26 




2.032 


3.822 


* 

300 




-64.47 


-62.60 




.173 


.341 


200 




-62.83 


-61.60 




.526 


2.294 


150 




-63.73 


-62.80 




.911 


2.889 


100 




-67.93 


-67.38 




1.496 


2.581 


70 




-71.50 


-70.29 




2.921 


3.026 


Number 














in 














class 




3 


10 









TABLE 13. Temperature means at mandatory levels at Thule, Greenland 
(winter), corresponding to temperature extremes of T (300) at the (b) 1-5^ 
class, and (c) S-IO^S class of probability. The population of class (a) was 
zero. The standard deviations and of observed temperatures in the (b) 
and (c) classes are shown in the right side of the table. 



26 



Ht. (km) 




FIG. 3. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Thule, Greenland, corresponding 
to the ly^ cold extreme set of forcing temperatures Tj(300) of Table 
13. The thin solid line depicts the January mean Thule sounding 
(1967-70). The dashed curves are identical to those of Figure 1. 
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Level 

(mb) 


Temperature-means (^C) 
in classes 


Std. deviations 
relative to class 


(°C) 

i-mean 




(a) 

^.01 


(b) 

T -T 

\0l .05 


(c) 

T -T 

.05 .10 


0 

a 


^b 


a 

c 


1000 


-9.87 


-12.17 


-14.51 


4.302 


3.849 


8.440 


850 


-5.30 


-15.31 


-17.09 


2.269 


3.179 


5.024 


700 


-12.67 


-19.73 


-23.68 


0.806 


1.608 


4.870 


500 


-29.97 


-34.99 


-38.21 


1.304 


2.732 


3.773 


300 


-57.03 


-59.19 


-60.40 


0.695 


1.041 


2.115 


★ 

200 


-74.70 


-68.61 


-64.80 


.084 


1.607 


.480 


150 


-73.57 


-66.81 


-63.74 


0.173 


2.496 


2.216 


100 


-73.43 


-71.30 


-68.29 


0.733 


2.413 


2.339 


70 


-77.50 


-75.46 


-71.52 


1.145 


3.156 


2.767 


Number 

in 

class 


3 


7 


8 









TABLE 14. Temperature means at mandatory levels at Thule (winter), corres- 
ponding to temperature extremes of 1/200) at the (a) 1^ class, (b) 1-5^ 
class and (c) 5-lOf^ class of probability. The standard deviations O , 
of the observed temperatures within each class are shown in the right 
half of the table. 
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20 




FIG. 4. The heavy solid line shows the mean regress ion- determined 
vertical temperature sounding over Thule, Greenland, corresponding 
to the Vja cold extreme occurrences of T (200) of Table 14, The thin 
solid line depicts the January mean Thule sounding (1967-70). The 
dashed lines are identical to those of Figure 1. 
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Level 


Temperature-means (^C) 


Std. 


deviations 


(°C) 


(mb) 


in classes 


relative to class-mean 






(b) 


(c) 












T -T 

^.01 .05 


T -T 

.05 .10 






a 

c 


1000 




0.05 


1.48 




.350 


1.147 


850 




-6.65 


-4.07 




.650 


2.086 


700 




-16.15 


-13.09 




.250 


1.953 


500 




-31.30 


-27.79 




2.000 


4.086 


300 




-40.05 


-45.54 




2.250 


3.443 


* 

200 




-36.70 


-38.31 




.600 


.330 


150 




-38.90 


-39.67 




1.000 


.575 


100 




-40.95 


-40.47 




.350 


.860 


70 




-40.25 


-39.56 




1.150 


1.498 


Number 














in 














class 




2 


16 









TABLE 15. Temperature means at mandatory levels at Alert (summer), corres- 
ponding to temperature extremes of *^200) at the (b) 1-55^ class and (c) 5- 
10^ class of probability. The population of class (a) was zero. The 
standard deviations- CJ^, (J^ of the observed temperatures within classes h 
and £ are shown on the right side of the table. 
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FIG. 5. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Alert, NWT, corresponding to the 
1/^ warm extreme occurrences of T_(200) of Table 15. The thin solid 
line depicts the July mean Alert sounding (1967-70). The dashed 
curves are identical to those of Figure 1. 
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Level 

(mb) 


Tempera ture~mea ns (^C) 
in classes 


Std. deviations 
relative to class 


(°C) 

:-raean 




(a) 

T 


(b) 
T -T 


(c) 

T -T ^ 


0 


a. 


a 




.01 


.01 ,05 


.05 .10 


a 


b 


c 


1000 


0.00 


2.68 


2.29 


.030 


1.384 


2.734 


850 


-4.65 


-3.93 


-3.21 


2.650 


2.318 


2.067 


700 


-13.75 


-12.48 


-11.94 


2.650 


2.467 


2.689 


500 


-26.25 


-26.42 


-26.48 


3.050 


1.117 


2.317 


300 


-43.50 


-49.23 


-46.78 


1.200 


1.168 


3.177 


200 


-38.10 


-38.83 


-39.09 


.800 


.619 


.909 


★ 

150 


-38.20 


-38.80 


-39.57 


.300 


.224 


.226 


100 


-40.60 


-40.33 


-40.33 


.000 


.500 


.469 


70 


-39.45 


-39.93 


-39.43 


.349 


.472 


1.471 


Number 

in 

class 


2 


6 


14 









TABLE 16, Temperature means at mandatory levels at Alert (summer), corres- 
ponding to temperature extremes of Tj(150) at the (a) 1^ class, (b) 1-5^ 
class and (c) 5-105^ class of probability. The standard deviations O , 

O of the observed temperatures within each class are shown in the 
right side of the table. 
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FIG. 6. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Alert, NWT, corresponding to the 
Vjo warm extreme occurrences of Tj(150) of Table 16 o The thin solid 
line depicts the July mean Alert “^sounding (1967-70). The dashed 
curves are identical to those of Figure 1. 
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1 

Level 
( mb 


Tomperature-mcans (^C) 
in classes 


Std. deviations 
relative to class 


(°C) 

-caean 






(b) 

T -T 

.01 .05 


(c) 

T -T 

.05 .10 




^b 


a 

c 


1000 




0.90 


2.56 




.100 


2.821 


850 




-2.42 


-3.72 




3.544 


1.943 


700 




-10.10 


-11.43 




3.412 


3.283 


500 




-25.80 


-26.69 




4.513 


2.884 


300 




-45.25 


-46.81 




2.296 


2.772 


200 




-39.80 


-39.85 




1.404 


1.211 


150 




-39.80 


-40.86 




1.931 


.875 


* 

100 




-40.20 


-41.16 




.122 


.105 


70 




-39.92 


-40.82 




1.803 


.699 


Number 

in 

class 


s 

i 

1 


4 


14 









TABLE 17 . Temperature means at mandatory levels at Thule (summer) , corres- 
ponding to temperature extremes of Tj(100) at the (b) 1-5^ class and (c) 
5-10> class of probability. The population of class (a) was zero. The 
standard deviations of the observed temperatures within classes ^ 

and £ are shown on the right side of the table. 
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FIG. 7. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Thule, Greenland, corresponding 
to the 1/^ warm extreme occurrences of Tj(lOO) of Table 17. The thin 
solid line depicts the July mean Thule sounding (1967-70). The dashed 
curves are identical to those of Figure 1. 
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Table 15 «hov8 that class (b) warm atmosphere is warmer than the 
class (c) case in the layer 400 mb to 125 mb. Above 125 mb Fig. 5 shows 
that the deviation between the mean case and the listed extreme at any 
level is small. For example, the standard deviation at Alert at the 
100 mb level is 1.26^C (and is smaller still at 70 mb). Below 400 mb, the 
temperature of the class (b) 'extreme^^ atmosphere is consistently colder 
than that of the class (b) atmosphere^ All of these features 
are also to be observed in Fig. 15, which gives a comparison of the 1^ 
warm extreme atmosphere and the July mean (1967-70) over Alert. 

A similar set of conclusions applies to the regression-dependent 
class (a) atmospheres over Alert based upon = T(150) as forcing variable. 
Here the results included in Table 16 indicate that the 500 mb tempera- 
ture is a level of small Interdiumal variability, and is in fact the 
location of a crossover between the T(extreme) and T(mean) curves. Below 
500 mb, the inferred temperature structure corresponding to a forcing 
temperature at the 1^ warm extreme at p == 150 mb shows an atmosphere 
slightly colder than the July mean. 

The summer extreme atmosphere at Thule, based upon a forcing- level 
warm extreme in Tj “ T(IOO) is listed in Table 17 and also graphically 
depicted in Fig. 7. The results summarized in column (b) of Table 17 
show that the extreme warm case at 100 mb is associated with mandatory 
levels which are consistently warmer at all levels (other than the surface) 
than in the class (c) atmosphere. 

5 . Recommendat ions 

For the best estimate of the forced- level extreme atmospheres con- 
sidered here, one merely reads off the temperature of the left-hand extreme 
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class of Tables 11,..,, 17. To be more realistic, one should allow a tem- 
perature-range of one standard-deviation ± O about any class (a) extreme 

a 

atmosphere. However, if a single temperature-estimate must be used at each 

level, the mean T -profile as a function of P is to be recommended. If 
a 

class (a) is lacking, then the T, -profile should be used. 

b 

The work presented here does not include certain other levels listed 
in the preliminary MIL STD 210-B atmosphere. The four years of climatologi- 
cal data for the years 1967-70 for these additional stations (see Table 2) 
is still undergoing data analysis. The results for OJMJAKON, Siberia, have 
been analyzed statistically and have been found to determine regression- 
generated mandatory- level temperatures much like those derived for the 
North American stations. It is expected that a study of the additional 
stations listed in Table 2 will yield fruitful results. 
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APPETTOIX 



It is not feasible due to the limited sample sizes to employ a t- 
test on the difference of means for the 1^ class relative to the 1-5^ 
class of Tables 11,..., 17. However, a well-defined test of signifi- 
cance of class-means is afforded by the t-test applied to the differ- 
ence of sample-means obtained by the results of the predictions of 
the means of the 10^ extreme relative to the full sample predictions 
(e.g., the means contained in Tables 3, 4,..., 9). 

The t-test on the predictions of the full-sample and the 10^ ex- 
treme involves at every mandatory level and station a test for signi- 
ficance of the pooled t-statistic after prediction. 



tp(Hf I. 




(N + n-11} 



(n- 2 ) 0 ^ .10 



-,1/2 

] 



1/2 



(A-1) 



Here the T's are the sample temperature means at the mandatory levels 
based upon the identical prediction equations. These means are listed 
in columns 2 and 7 of Tables 3, 4,..., 9. C is the standard error of 
estimate listed in column 5, and 0_ that of coluoa *10, associated 

^ • .L VJ 

with the prediction of the extreme sample. The number of degrees of 
freedom associated with these G^-values are N-9 and n-2 respectively, 
where 



N = full-sample size 
n = nominal 10^ extreme sample size 

The pooled tp-parameter is then tested using ^ + n-i;^ degrees of freedom 
for tp. 
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At the 70-mb level in Table 3, the tp-statistic based upon the 
listed difference of the means is tp = 3.113. Using the number of 
degrees of freedom N + n-11 = 258, the sample tp-statistic is signi- 
ficant at a level of 99 . 9 % probability. 

In a similar manner all sets of differences of full-sample and 
"extreme" sample means obtained by the prediction method of this 
study prove to be significant at levels in excess of 99 % confidence. 
This may be verified by applying all the statistics of Tables 3, A, 
...,9(a,b) to the test-statistic tp o£ Eq. (A-1). 

In this Appendix, t-test procedures for considering definitive 
tests of significant differences of means between full- and extreme- 
sample statistics at mandatory levels have been set forth. In 
defining the extreme-sample, the 10^ probability extreme at the 
forcing level has been considered as the basis of the sampl^ includ- 
ing the resulting regression-statistics. Beyond the 10^ probability 
level, reduced sample sizes preclude the usual t-test procedures. 
However, the stratification 5 based upon extreme sub-samples (a), (b) , 
(c) of Tables 11, 12,..., 17 give rise to regression-generated 

T (meai^ which might well have been extrapolated from the 10^ sub- 

Si 

sample results of Tables 3, 4,..., 9, respectively. 
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